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We apply a recently developed stochastic method to the Shastry-Sutherland model on 
4x4 and 8x8 lattices. This method, which we call the Stochastic State Selection Method 
here, enables us to evaluate expectation values of powers of the Hamiltonian with very 
limited number of states. In this paper we first apply it to the 4x4 Shastry-Sutherland 
system, where one can easily obtain the exact ground state, in order to demonstrate that 
the method works well for this frustrated system. We numerically show that errors of 
the evaluations depend much on representations of the states and that the restructured 
representation is better than the normal one for this model. Then we study the 8x8 
system to estimate energy eigenvalues of the lowest S = 1 state as well as of the lowest 
excited 5 = state, where S denotes the total spin of the system. The results, which 
are in good accordance with our previous data obtained by the Operator Variational 
method, support that an intermediate spin-gap phase exists between the singlet dimer 
phase and the magnetically ordered phase. Estimates of the critical coupling and of 
the spin gap for the transition from the dimer phase to the intermediate phase are also 
presented. 
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1. Introduction 

Considerable effort has been exerted to develop effective methods to calculate various quantities 
in quantum spin systems. 1-3 ) A new method we have proposed in a previous paper 4 ) is a kind of 
Monte Carlo approach, where stochastic variables play an important role. It is completely different 
from the conventional quantum Monte Carlo methods which employ random walks or importance 
samplings. 3 ^ In our method, which we call the Stochastic State Selection (SSS) method hereafter, 
random variables are used to reduce the number of states in the vector space which is huge for 
most systems of large sizes. To generate these random variables we introduce a probability function 
named on- off probability function, which "switches off 1 a number of states in the vector space so 
that we can calculate approximate expectation values of powers of the Hamiltonian from a small 
number of the "on" states. Repeating this process with a new set of random variables we can 
obtain averaged expectation values which are very close to the exact ones. 

The fundamental idea in the SSS method is deeply connected with the variational approach 
to quantum spin systems, where one searches for an effective basis, a set of states which can be 
treated within the limit of computer facilities. In the variational approach there are two essential 
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steps: increasing the number of states in the wave function under current investigation, and then 
reducing it. In the former step one usually operates the Hamiltonian to the wave function, which 
is a mathematically established procedure. How to formulate the latter step, on the other hand, is 
quite controversial. In early works 5 ) states with small coefficients are dropped, but the resultant 
wave functions are not satisfactory. Sampling the states using the Monte Carlo techniques has 
also been suggested in recent works. 6 ' 7 ) In the SSS method one can reduce the number of states 
based on a mathematically justified way. It is this remarkable property that makes us to believe 
the method is worth investigating. 

In the previous paper we described the SSS method and applied it to the two-dimensional spin 
one-half Heisenberg model. In this paper we demonstrate that it can be applied to the system 
suffering from the negative sign problem, too. A concrete example is the Shastry-Sutherland (SS) 
model 8-17 ) on a 4 x 4 lattice, whose results reproduce the exact values in good precision. Next we 
try to obtain physical quantities, the energy eigenvalues with fixed values of total spin S, for the 
8x8 lattice changing the coupling ratio across the critical region. 

In the next section we briefly explain the SSS method. After defining the on-off probability 
function and the random choice matrix, we describe how we calculate expectation values of H L , H 
being the Hamiltonian of the system under study. Section 3 is a short description of the SS model, 
where we explain the Hamiltonian and several properties of this model. In section 4 numerical 
results on a 4 x 4 lattice near the critical point are presented for two representations, normal 
and restructured. 1 ^ Deviations of the data here indicate that the restructured representation is 
more effective for the coupling ratio we study in this section, which supports the picture that 
the orthogonal dimers dominate here. Based on these discussions we try in section 5 to obtain 
physical quantities on a 8 x 8 lattice, the energy eigenvalues with the total spin S = and 5 = 1, 
in the intermediate region of the coupling ratio. As noted previously, 4 ) there are two model- 
dependent problems to be solved. One of them is how we generate a good approximate state | i/ja)- 
Another is how we extract eigenvalues from the expectation values {ipx I H l \ ipx). In addition 
to the calculations of expectation values for the restructured representation, extensive discussions 
on these subjects are given in this section. We will show that we successfully obtain the results 
on the energy eigenvalues which indicate a phase transition from the singlet dimer phase to an 
intermediate spin-gap phase. Finally, the last section is devoted to summary and discussions. 

2. Stochastic State Selection Method 

Here we give a brief description of the SSS method. 4 ) We define the on-off probability function 

by 

P(v) = ls( V -a) + (l-^)S(r,), (1) 
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where a is a constant which is greater than or equal to 1. It is clear by definition that 

/>oo 

((1)) EE / P( V )d V = l, (2) 
J 

«r/» ee / r,P(r,)d V = 1, (3) 
J 

((ri 2 )) = / rj 2 P(r,)dr, = a, (4) 
J o 

where (()) denotes the statistical average. 

Let us expand a state | ip) by a basis { | i}} 



N v 

i=i 



(5) 



where Ny denotes the size of the full vector space, and define the expectation values E(L) (L 
1,2, - ••) as 



E(L) = (ip \ H L \ ip) . 
We then introduce the random choice matrix, an Ny x Ny diagonal matrix, 

( 771 ••• \ 

o m ■■■ o 



(6) 
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■■■ VN V / 
Each random variable rji in (7) is determined by the on-off probability function 



V 



(7) 



Pi(m) = —Km - a i) + (! - — )<5(^), 



where 



a,- = 



max(l,e/|ci|) (if q / 0) 
e/<5 (if Q = 0) , 



(8) 



(9) 



with given constants e and 8 (0 < <5 < e). Using (L + 1) independent random choice matrices 



M {r)(m)} = diag.{ri { ™\r] { ™\- ■ ■ ,r?^} 19) we define 



E {v} (L) ee (V I M WL+1)} HM {v(L)} HM {v(L - 1)} ■ --HM^y \ 



(10) 



It is easy to see that ((E{ v y(L))) = E(L). Note that by operating M^ m )y we can obtain a truncated 
vector space for HM^ m -i)y • • • HM^^y \ ip), since r\i^ is zero with the probability 1 — 1/dj. Also 
note that the variances of rlf^ are the same for all m because both ((rj^)) and (([r]^} 2 )} do not 
depend on m. 
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In numerical studies we measure 

^smpl 



«£ W (L)» smp i = E {v}k (L), (11) 



k=l 

p\ v} {L) = (([E {v} (L)] 2 )) smpl - ((E {v} (L))f smpl , (12) 
using n smp i samples. The error of (L) is evaluated by 



£ r(L ) = 2 ^k^. (13) 

V ^smpl 

The fact that the number of non-zero components much reduces by operating M^ m )^ to 
HM^( m -i)y ■ ■ ■ HMj^i)-! | ip) is essential to numerical calculations. Let us denote by ((-/V&(m))) smp i 
and {{N a (m))) sm pi the number of non-zero coefficients we count in the measurement before and 
after we operate M^( m )y to HM^ m -i)y ■ ■ ■ HM^ m y \ ip), respectively. 

3. Shastry-Sutherland Model 

As is discussed in several papers, 9-17 ) one of recent fascinating topics in two-dimensional quan- 
tum spin systems is the Shastry-Sutherland (SS) model, 8 ) the model of orthogonal dimers with the 
intra-dimer coupling J (> 0) and the inter-dimer coupling J' (> 0). 

In this model a spin is located at the sites 

d d d d 

(2n x a+ -,2n y a + -), (2n x a - -,2n y a - -), 

{{2n x + l)o + ^, (2n y + l)a - |), ((2n x + l)a - ^, (2n y + l)a + ^), 

on the 2N x 2N lattice in Fig. 1, where n x = 0, 1, • ■ ■ , iV — 1 and n y = 0, 1, • • • , N — 1, 2a is the 
unit distance between dimers and \[2d is the distance of spins in a dimer. The Hamiltonian is given 
by 

N-l 

4 J 



1 N ~ 1 

H = -J ^2 {h a (n x ,n y ) + h b (n x ,n y )} 



n x ,n y =0 
N-l 



+ jJ' ^2 {hi(n x ,n y ) + h 2 {n x ,n y ) + h 3 (n x ,n y ) + h 4: (n x ,n y )}, 

n x ,n y =0 

where the partial Hamiltonians are, as given in ref.17, 



d d d <i, 

h a {n x ,n y ) = a(2n x a + -,2n y a + -) ■ cr(2n x a - -,2n y a - - y 



and so on with a(x,y) denoting the Pauli matrix at the location (x,y). 

This model provides a nontrivial system which relates an exactly solvable spin-gap model 
obtained in the J' — ► limit to a square lattice model where intra-dimer interaction vanishes 
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(J = 0). It is known that the exact ground state of the SS model is the direct product of the 
singlet dimers when J' < J' c ~ 0.68 J, 12 ) while for sufficiently large J' /J the system is in the 
magnetic order phase with the Neel ground state. It is almost doubtless that another phase exists 
between these two phases, but the location and the properties of this phase are not definite yet. n ~ 17 ) 
This phase is suggested to be gapless or nearly so in refs.13 and 15, while authers in refs.ll, 12, 
14, 16 and 17 conclude it is gapped. Numerical data for the spin-gap is presented only in refs.12 
and 17. 

In study of these phases a fundamental quantity to be examined is the difference between the 
energy of the singlet-dimer state Eq and the lowest energy £ as a function of the coupling J'/ J, 
when the momentum p=(0, 0), the total spin S = j (j = 0, 1) and the /+ parity is even or odd, 
where the I + parity is the parity for the x-y reflection I + cr(x, y)l\. = cr(y,x). 

4. Numerical results on a 4x4 lattice 

In this section we show the numerical results obtained on a small lattice in order to demonstrate 
that we can apply the SSS method to the SS model with J' ~ J' c , where the negative sign problem 
is serious, without any difficulty 

It is quite easy to calculate the exact ground state | ips) of the model on the 4x4 lattice using 
the Lanczos techniques. We therefore do not need to pursue any trial state here. We just measure 

E E { V } (L) = (Ve I M^L+vyHM^wyHM^L-v y HM {vW } | Ve) (14) 

(L = 1, 2, • • • ) and compare them with Ee(L) = E L , where E denotes the eigenvalue of H, namely 
H | Ve> =| ife)E. 

In this section we fix the coupling ratio J'/J = 0.7, a value slightly above the critical point, 
for which E/J = —6.3502. Here we calculate ((E E { v y(L))} smp i using two representations. One 
of them is the normal one, where the state on each site of the lattice is represented by the z 
component of the spin sitting on the site. In another representation we named restructured, two 
spins connected with the coupling J are rearranged in three triplets and one singlet states. 18 ^ 
The total number of non-zero coefficients in expansion of | Ve) is 12,870 (8,565) in the normal 
(restructured) representation. 

Table I shows our results up to L = 10 obtained from 10 4 samples with e = 0.1, together with 
the exact values. Note that 5 in (9) is not necessary for | Y>e)- We see that the results using the 
restructured representation are in good agreement with the exact ones. Especially for small values 
of L the agreement is excellent. The relative errors are in the range from 0.11% (L = 1) to 2.8% 
(L = 10). The results from the normal representation are also compatible with the exact values. 
We, however, observe that the statistical errors are much larger than those in the restructured 
representation except for very low L. The relative error is 0.12% for L = 1 but it rapidly increases 
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as L grows to reach 113% (L = 9), where the evaluated value is not meaningful. For larger values of 
e such tendency is more outstanding. With e = 0.2 and n smp i = 10 4 , for example, we cannot obtain 
reliable values of ((£ , E{ ?? }(i)))smpi (L > 6) in the normal representation, while the results from the 
restructured representation are satisfying up to L = 10. When the number of samples is increased 
we are in principle to obtain satisfactory evaluations even in the normal representation. In fact, 
we observe that for 10 5 samples with e = 0.1 the statistical errors become smaller than those for 
10 4 samples except for L = 10. ((E E ^(9))) smp i/ J 9 and ((£ ; E { J? }(10))) smp i/J 10 for 10 5 samples are 
-(0.189 ± 0.066) x 10 8 and -(0.120 ± 0.088) x 10 9 , respectively. 

Table II presents the numbers of non-zero coefficients measured in the same calculations with 
10 4 samples, which is to show how the number of non-zero coefficients decreases by each operation 
of the random choice matrix. In both representations we see that ((N^(L)}} smp \, which denotes the 
number of the "on" states after the operation, is less than a tenth part of ((N^(L))) smp \, the number 
of non-zero coefficient before the operation. In addition, the data suggest that both ((N^(L))) smp i 
and ((Nf(L))) smp \ become constants for large L in each representation and these constants very 
weakly depend upon the representation. 

Results in Table I strongly indicate that the restructuring is quite effective here. What is the 
reason for such qualitatively different behaviors in these two representations then? To answer this 
question we plot in Fig. 2 the distributions of sorted coefficients in the expansion of | Ve) hi both 
representations. The abscissa is the basis number i' which is reordered so that |cj/| > \cf \ holds 
for i' < f. We observe the distributions are different in the region of the small coefficients, where 
the normal representation has a longer tail. The minimum value of the non-zero coefficients in the 
normal representation is 7.6 x 10 -6 while it is 1.5 x 10~ 4 in the restructured representation. We 
will return to this issue in an Appendix. 

5. Numerical results on a 8x8 lattice 

Let us proceed to a larger lattice with 8x8 sites and show that the SSS method is helpful to 
evaluate the energy of the first excited states. We study the system with the total spin 5 = and 
1 in the parameter region 0.55 < J'/ J < 0.8. Here we use 

E {v} {L) = <V> | HM {v(L)} HM {v(L -i )} • • • HM {vW} | V) (15) 

instead of (10) in order to decrease the statistical errors. Throughout this section we employ 
the restructured dimer representation, which is a very powerful technique for the model under 
investigation. The results we present here are compared to those in ref.17, where we have been 
successful to obtain approximate values on the energy difference for this system by means of the 
operator variational (OV) method. 

For this lattice size it is necessary to find an approximate state | ipA) without knowing the 
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exact eigenstate | Ve)- Here we employ a systematic way using the Lanczos-like method in small 
vector spaces; we start with a trial state | Vtri), calculate the matrix elements of H from the basis 
{ | Vtrl),# | Vtrl),-" I Vtrl),--- > H p | Vtrl)} with a small integer p, numerically diagonalize the 
matrix to obtain its all energy eigenvalues and define | Va) as the eigenstate which provides the 
lowest eigenvalue among them. In the following calculations we fix p = 6. 

As for the trial state | V'trl) we employ one of the states described in the appendix D in ref.17 
whose energy eigenvalue is found to be the lowest of all states we tried in the numerical study by 
means of the OV method. Namely we start from | Vtrl) = I *nn.,o,(o,o),i)') for S" = and | V'trl) = 
I ^nn,i,(o,o),d') f° r S = 1, which are among the two-nearest-neighboring-triplet states with the odd 
and the even I + parity, respectively. It should be noted that in our previous study 17 ) the energy 
eigenvalue for the even I + is much higher than that for the odd I + when S = while for S = 1 the 
lowest energy is independent of the I + parity. 

Now we present the results. Table III and IV show values of ((E^ v j(L))) smp i observed for 
the S = and S = 1 cases, where the suffix A represents that | V) in (10) is an approximate state 
| Va)- In both tables we see very small errors, which support that the restructuring technique is 
quite effective. How many non-zero coefficients appear here? When S = the number of non- 
zero coefficients in | Va) is 1,264,256, while ((N^(5))) smp \, the number of non-zero coefficients in 
i^M| r) (4)}i^M| J) (3)}i^M| J) (2)}^M| ?? (i)| | Va) in measurements, ranges from (9.78 ± 0.03) x 10 5 at 
J' I J = 0.55 to (1.912 ± 0.004) x 10 6 at J'/J = 0.75. In the S = 1 case | Va) contains 690,624 
components and we observe ((iV 6 A (5))) smpl is (2.84±0.02) x 10 5 ((1.534±0.004) x 10 6 ) at J'/J = 0.55 
(0.8). By operating M^^y we can reduce it to less than a tenth. In concrete terms, ((-/V" A (5))) smp i 
for S = ranges from (8.35 ± 0.02) x 10 4 at J'/J = 0.55 to (1.657 ± 0.003) x 10 5 at J'/J = 0.75 
and for S = 1 it is (2.27±0.01) x 10 4 ((1.271 ±0.003) x 10 5 ) at J'/J = 0.55 (0.8). Thus we see that 
we can evaluate expectation values of H L not only for small J'/ J but also in the critical region of 
the parameter space. 

How do we extract energy eigenvalues from these data? We try two ways, one is to use the 
extrapolation formula and another is to do the Lanczos-like calculation with the basis { | Va), 
H | Va), H 2 I Va)}, noting that the Hamiltonian matrix elements from this basis are given by 
(Va I H L | Va) up to L = 5. 

The extrapolation formula is the one discussed in the previous paper, 4 ) 

(Va I H L | Va) = E L (q + J^ 1 - - ) = F(L, E, q , gl , a) . (16) 

L + a + 1 

In process of determining four parameters in this formula, E, qo, q\ and a, we use two relations 

F(0,E,q ,qi,a) = q + — ^— = (Va | Va) = 1 , 

a + 1 

F(l,E,q , qi ,a) = E ■ (q + ) = (Va | H \ Va) , 
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so that qi and a are calculated from E and qo- Note that the value of (Va | -H" | ^a) is known in 
the process of generating | Va)- Then we look for the values of E and qo which minimize 

2 



L=l 



1 _ ((- E A{r ? }(-^)))smpl 

F(L,E,q ,qi,a) 



(17) 



within the ranges (Va | H \ Va) — 0.5 < E < (Va \ H | Va) and 0.5 < qo < 1.0. Results of the 
fitting for J'/ J > 0.5 are shown in Tables V (S = 0) and VI (5 = 1). The statistical errors of the 
measured values, 



52Er = 



5 

E 

L=l 



Er(L) 



1 2 



(( E A{r,}(L)))smpl 



(18) 



are also presented in the tables. The fact that S&t is smaller than S^Er means the fitting is 
statistically acceptable. We see that values of qo obtained from the fit are close to 1. Since qo 
represents the overlapping between the trial state and the exact state, (Va I Ve) ; this implies that 
our way to generate the state | Va) supplies good approximations of the exact eigenstate | Ve) for 
this model. The energy difference from the fit normalized by the intra-dimer coupling, (E{^ — Eq) / J, 
is plotted in Fig. 3 (S = 0) and Fig. 4 (S = 1) together with those obtained by the OV method. 
Here Eo denotes the energy of the ground state in the singlet dimer phase, which is —24.0 J for the 
8x8 system. We see that the data from the OV method with L max = 4 and the data from the SSS 
method are in good agreement when J' /J = 0.55. For larger coupling ratios the data from the SSS 
method are below the data from the OV method. Especially near the critical point this tendency is 
outstanding, suggesting that the SSS method is more suitable to study the phase transition in the 
SS model. The data by the SSS method also support that an intermediate spin-gap phase exists, 
because at J'/ J ~ 0.71 E^ t — Eo becomes zero for S = while E& t — Eo ~ 0.18 J for S = 1 there. 

Let us turn to another way we try to extract energy eigenvalues. Based on the observation that 
| iPa) is a good approximate of | ^e), we apply the Lanczos-like method to evaluate the eigenvalue 
j^( n ) ( n= i 5 2, 3) from the matrix whose matrix elements are calculated from bases {(V'a)}) 
{| ip A ) ,H I ^ A )} and {| tp A ) ,H \ V'a) , H 2 \ ^a)}, respectively. Substituting ((^A{r,}(^)))sm P i in 
Tables III and IV for E A (L) = (^a | H l \ Y> A ) in 

ai = (V'A | H | Va) , 



Pi = \l (V'A | H 2 | Va) - a? 



= -m ((V'A I H 3 | Va) - 4) ~ 2ai , 



° 2 - 01 



P2 = 



^Jj- ((Va I H 4 | Va) - 2«i(Va I # 3 I Va) + «?(Va I H 2 \ Va)) -a 2 - 0\ 
"3 = i (W>a I H 5 | Va) - 2«i(Va I H 4 | Va) + a?(V A I H 3 \ Va)) - 2a 2 
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-| («i + 2q 2 ) 

P2 



«2 
J, 



2 ' 



we diagonalize 






/ 


Oil 




\ 








«2 






V 







«3 / 



tf« = (ai), H™ = 



The results are also in Tables V (S = 0) and VI (S 1 = 1). As we see in tables we obtain better 
estimates by E^ than by Note that values of E^ are in good accordance with E^ t in 

Tables III and IV, which indicates that the fits are reliable. We could not, on the other hand, 
obtain trustworthy values of E^ from the present data. It is mainly because of the insufficient 
statistical precision, which we can cope with by increasing the number of samples, or by decreasing 
the parameter e. It seems that the latter is much more effective than the former. In the 5 = and 
J'/ J = 0.7 case, for instance, we observe that we cannot obtain any E^> from 9 x 10 4 samples with 
e = 0.001 while E^ calculated from 10 4 samples with e = 0.0005 is -23.928J, which is slightly 
lower than the value of E^ in Table V. 



6. Summary and Discussions 

In this paper we calculated the expectation values of H L (L = 1, 2, • • • ) of the two-dimensional 
SS model 8 ) using the SSS method. 4 ) Comparing our results on a 4 x 4 lattice with exact values 
we demonstrated that the method is applicable even to this strongly frustrated quantum spin 
system. We also showed that deviations of the evaluated values are much less in the restructured 
representation than in the normal representation. In the normal representation of the SS model 
we observed rapid increase of the statistical errors for large values of L, which is in contrast with 
our previous study of the 4x4 spin one-half Heisenberg model. A > A qualitative discussion on the 
reason why such increase occurs is given in the appendix. 

In study of the 8 x 8 SS system we concentrated our attention on an intermediate phase reported 
in refs. 11-17. We examined energy eigenvalues of the lowest excited states with the total spin S = 
and S = 1 for several values of the coupling ratio. It should be emphasized that the method is 
powerful even in the critical region of the model, where the negative sign problem prevents us 
from obtaining statistically meaningful results by means of the standard quantum Monte Carlo 
techniques. 20 -* In order to numerically extract the energy eigenvalues from the observed {H L ) we 
tried two ways, a fitting within some parameter space and a Lanczos-like diagonalization. We see 
results from these two ways are consistent to give us reliable estimate of energy eigenvalues. It 
should be kept in mind that employing highly qualified trial states in the SSS method is a key for 
the success. Thus we see that the critical value of J' / J is less than 0.71 for the phase transition 
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from the singlet dimer phase to the intermediate phase, which is compatible with the value in 
ref.12. Note that this upper bound is obtained with a non-perturbative method. The spin gap at 
the critical point, ~ 0.18J, is about one-half of the one reported in ref.12, however. More intensive 
investigations should be made in future work before we establish physical properties of this phase. 

Appendix: Large Variances 

In this appendix we discuss the reason why very large variances are observed in the normal 
representation of the SS model. As we have seen in Fig. 2 much more terms with very small coeffi- 
cients appear in the expansion of the wave function | in the normal representation compared to 
the restructured representation. So we analyze the expression for the variance paying our attention 
to contributions from those terms. 

Let us examine <r^(L), the variance of ((E{ v y(L))) introduced in ref. 4, 



a\ n} {L) ee ((E {n} {L) 2 )) - ((E {V} {L))) 



■ hk'i'Cv 



= EE^'EE ^'"'^^ 

= E E ■■ ■ E E Cihij ■ ■ ■ hkici ■ Ci'hi'j' ■ ■ ■ h vv c v ■ {1 + 8 iV {{{r} 2 )) - 1)} 



3,3 



x {1 + -!)}•••{! + S kk/ (((vD) ~ !)}{! + M«u?» " 1)} 

E(Lf . 



(Al) 



For L > 2 this variance contains one or more terms that might blow up with small non-zero 
coefficients, because the statistical average ((r/ 2 )) = e/|c s |, which follows from (4) and (9), becomes 
huge if a non-zero coefficient c s in (8) is much less than e . The most dangerous term among them 
is 



T = J2 &l(hi j ?(h ]k ) 2 ---(h lm ) 2 {h mn ) 2 

i,j,k,— ,l,m,n 

x «k 2 » - l)(«f£» - l)(«u2» -!)••• (((rif)) ~ 



!)(«)- I)- 



For L = 1 this term is harmless since in this case 



^2 c i m & x ( 1) i — r ) — 1 ^2 c i max ( 1) I — r 

c^o L V l c «l / J c .^lq V \ \ c 3\ 

E ( e M _c * 2 ) E ( e l c il ~ c i)( /l u) 2 • 

0<|e;|<e 0<|cj|<e 
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Actually we observe on a 4 x 4 lattice that the measured variances of ((E{ v y(l))), which are defined 
by (12), are comparable in both representations. They are 0.153J 2 and 0.135J 2 when e = 0.1 and 
^smpi = 10 4 for the normal representation and the restructured representation, respectively. 
Let us then consider the case L = 2, where 

i k j 

= E ( £ M" c2 ) E (^k\-4)J2(hij)\h jk )\((r,]))-l) 

0<|C 1 |<£ 0<|c fc |<£ j 

~ e ( £ M- c2 ) E ^c k \-4) e (M 2 (M 2 A- 

0<|ci|<e 0<|c fc |<e 0<|cj|<e J ' 

If the Hamiltonian matrix is positive-definite (hj k > 0) and all coefficients of the eigen state are 
non- negative (c k > 0), 

< h jk c k < E hjm = Ecj 
i 

holds for any k. In this case it follows that 

T ~ E ( ec i- c h E ( e_c fc) E ( h ij) 2h jk- hjkCk' — 

0<Ci<e 0<c k <e 0<Cj<e 3 

< E ( ec i~4) E ( e ~ c k) E {hijfhjk ■ E Cj ■ ^~ 

0<Ci<e 0<c k <e 0<Cj<e 3 

= eE E ( ec i - c i) E ( e ~ Cfc ) E ( h ij) 2fl jk , 

0<Ci<e 0<c k <e 0<Cj<e 

which gives a safe upper bound of T. This means that the positive-definite systems are free from 
unacceptably large variances caused by the smallness of the coefficients. For the frustrated systems 
with the negative sign problem, on the contrary, there is no guarantee that they can escape from 
this difficulty. It is very likely, therefore, that the frustration and a lot of small coefficients would be 
responsible together to such behavior of the statistical errors in the measurement. Clearly similar 
discussions are possible for L > 3. 
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Normal 




Restructured 




L 


E E (L)/J L 


«£ EM (L)» 


/ T L 
smpl/ J 




1 


-0.63502 


xlO 1 


-(0.63484 ± 0.00078) 


xlO 1 - 


-(0.63486 ± 0.00073) 


xlO 1 


2 


0.40325 


xlO 2 


(0.40336 ± 0.00071) 


xlO 2 


(0.40298 ± 0.00063) 


xlO 2 


3 


-0.25607 


xlO 3 


-(0.25580 ± 0.00096) 


xlO 3 - 


-(0.25601 ± 0.00052) 


xlO 3 


4 


0.16261 


xlO 4 


(0.1619 ± 0.0023) 


xlO 4 


(0.16277 ± 0.00052) 


xlO 4 


5 


-0.10326 


xlO 5 


-(0.1056 ± 0.0052) 


xlO 5 - 


-(0.10307 ± 0.00053) 


xlO 5 


6 


0.65576 


xlO 5 


(0.692± 0.043) 


xlO 5 


(0.6568 ± 0.0042) 


xlO 5 


7 


-0.41642 


xlO 6 


-(0.425 ± 0.053) 


xlO 6 


-(0.4148 ± 0.0035) 


xlO 6 


8 


0.26443 


xlO 7 


(0.277 ± 0.081) 


xlO 7 


(0.2655 ± 0.0034) 


xlO 7 


9 


-0.16792 


xlO 8 


-(0.06 ± 0.19) 


xlO 8 


-(0.1701 ± 0.0025) 


xlO 8 


10 


0.10663 


xlO 9 


(0.107 ± 0.048) 


xlO 9 


(0.1048 ± 0.0030) 


xlO 9 



Table I. Results on ((S B {^}(L)}) smp i (L = 1, 2, • • • , 10) obtained for the 4x4 Shastry-Sutherland (SS) model by the 
Statistic State Selection (SSS) method. The number of samples is 10 4 and e = 0.1 for both of the normal and the 
restructured representations. Exact values of E-e(L) are also presented for comparison. 
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Normal Restructured 



T 

Li 


// atE( t\\\ 


// atE( t\\\ 


// atE/ t\\\ 


// /vE ( T\\\ 


1 


IZO ( U 


KTV 1 
1 1 . 1 


OOOO 


041. U 


O 

z 


( ZOU.Z 


400.0 


1 404. 1 


^QO 8 
0OZ.8 


q 
o 


COM V 


Ada n 

4DO.U 


04»4.U 


oUZ.o 


4 


oloy.U 


4DO.O 


^i en n 

oioy.u 


4yz.o 


c 



OilU.O 


40Z.4 


ou / y. / 


A 88 Q 

4oo.y 


6 


6101.0 


462.2 


6037.5 


487.5 


7 


6097.7 


462.3 


6020.0 


487.0 


8 


6099.1 


461.9 


6015.6 


486.6 


9 


6097.0 


462.0 


6010.9 


486.5 


10 


6097.7 


462.0 


6009.5 


486.3 



Table II. Numbers ol non-zero coefficients before and after operating the random choice matrix M^ L )^ to the 
state HM^L-i)} ■ ■ ■ HM^i)^ \ ips) for the 4 x 4 SS model obtained from 10 4 samples with e = 0.1. 
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L 










J / J = U.00 




J'/J = 0.60 




1 


— ^U.ZoZOOOOO ± U.UUUUUUZf^ 


xlO 2 


-(0.23467222 ± 0.00000030) 


X 1U 


2 


(0.5414619 ± 0.0000011) 


xlO 3 


(0.5507468 ± 0.0000012) 


xlO 3 


3 


-(0.12600061 ± 0.00000047) 


xlO 5 


-(0.12925980 ± 0.00000050) 


xlO 5 


4 


(0.2932180 ± 0.0000036) 


xlO 6 


(0.3033880 ± 0.0000032) 


xlO 6 


5 


-(0.682344 ± 0.000021) 


xlO 7 


-(0.712097 ± 0.000032) 


xlO 7 




J'/J = 0.65 




J'/J = 0.70 




1 


-(0.23680656 ± 0.00000033) 


xlO 2 


-(0.23908546 ± 0.00000036) 


xlO 2 


2 


(0.5608338 ± 0.0000014) 


xlO 3 


(0.5717199 ± 0.0000015) 


xlO 3 


3 


-(0.13283432 ± 0.00000055) 


xlO 5 


-(0.13673295 ± 0.00000061) 


xlO 5 


4 


(0.3146405 ± 0.0000031) 


xlO 6 


(0.3270462 ± 0.0000032) 


xlO 6 


5 


-(0.745321 ± 0.000023) 


xlO 7 


-(0.782349 ± 0.000030) 


xlO 7 


J'/J = 0.75 


1 


-(0.24149807 ± 0.00000039) 


xlO 2 






2 


(0.5833748 ± 0.0000016) 


xlO 3 






3 


-(0.14095261 ± 0.00000067) 


xlO 5 






4 


(0.3406212 ± 0.0000035) 


xlO 6 






5 


-(0.823245 ± 0.000028) 


xlO 7 







Table III. S = results on ((-B A {^}(L)» smp i (L = 
with n smp i = 10 4 , e = 0.001 and S = 0.00001. 



1, 2, • ■ ■ ,5) obtained for the 8 x 8 SS model by the SSS method 
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L 










J 1 J — I). 00 




J'/J = 0.60 




1 


/n OQ/inQd o _i_ n nnnnnm 
— t^U.Zo4Uoolo ± U.UUUUUUlo, 


xlO 2 


-(0.23504307 ± 0.00000016) 


\^ 1 n2 


z 


I^U.04/ /44/D zb U.UUUUUUOO, 


xlO 3 


(0.55246007 ± 0.00000065) 


X 1U 


Q 

o 


— (U.lZoiybUl ± U.UUUUUUzu, 


xlO 5 


-(0.12985518 ± 0.00000029) 


X 1U 


A 

4 


(U.oUUUoOO ± U.UUUUUZ1, 


xlO 6 


(0.3052268 ± 0.0000022) 


X 1U 


r 



( n 7fiooi o _i_ n AfifinoQ 1 * 
— (U./Uzzlo ± U.UUUUzo^ 


xlO 7 


-(0.717444 ± 0.000027) 


X 1U 




J 1 J — U.OO 




J'/J = 0.70 




1 
1 


fn oQfiooni/i -i- n nnnnnm ga 
— ^U.ZoOZZU14 zt u.uuuuuuiy^ 


xlO 2 


-(0.23762394 ± 0.00000024) 


X 1U 


2 


(0.55801950 ± 0.00000081) 


xlO 3 


(0.5647149 ± 0.0000010) 


xlO 3 


3 


-(0.13182375 ± 0.00000034) 


xlO 5 


-(0.13421698 ± 0.00000042) 


xlO 5 


4 


(0.3114219 ± 0.0000023) 


xlO 6 


(0.3190188 ± 0.0000023) 


xlO 6 


5 


-(0.735702 ± 0.000014) 


xlO 7 


-(0.758308 ± 0.000017) 


xlO 7 




J'/J = 0.75 




J'/J = 0.80 




1 


-(0.23936605 ± 0.00000029) 


xlO 2 


-(0.24150001 ± 0.00000036) 


xlO 2 


2 


(0.5731340 ± 0.0000012) 


xlO 3 


(0.5835473 ± 0.0000015) 


xlO 3 


3 


-(0.13726176 ± 0.00000049) 


xlO 5 


-(0.14106434 ± 0.00000057) 


xlO 5 


4 


(0.3287936 ± 0.0000028) 


xlO 6 


(0.3411166 ± 0.0000030) 


xlO 6 


5 


-(0.787688 ± 0.000021) 


xlO 7 


-(0.825086 ± 0.000022) 


xlO 7 



Table IV. 5 = 1 results on ({E AM (L))) smpl (L = 
with n smp i = 10 4 , e = 0.001 and 8 = 0.00001. 



1, 2, • ■ ■ ,5) obtained for the 8 x 8 SS model by the SSS method 
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J'lJ 




<S*2Er 




Sfit 


<Z0 fit 


E&/ J 


E^/J 


E&/J 


E^/J 


0.55 


1.1 


xlO" 9 


3.7 


xlO- 10 


0.995 


-23.276 


-23.269 


-23.272 


-23.273 


0.60 


2.1 


xl0~ 9 


8.4 


xlO" 11 


0.994 


-23.478 


-23.467 


-23.473 


-23.476 


0.65 


1.1 


xl0~ 9 


8.5 


xlO" 11 


0.988 


-23.700 


-23.681 


-23.691 


-58.851 


0.70 


1.6 


xl0~ 9 


1.7 


xlO- 10 


0.979 


-23.942 


-23.909 


-23.927 




0.75 


1.3 


xl0~ 9 


3.6 


xlO- 10 


0.975 


-24.196 


-24.150 


-24.178 


-24.183 


Table V. 


Results from the fit and from the Lanczos-like evaluation for S= 


=0 data in Table III. 


J'jJ 




S-mx 




S&t 


<7o fit 


E&t /J 


E^/J 


E^/J 


E^/J 


0.55 


1.1 


xlO' 9 


4.9 


xlO" 11 


0.986 


-23.410 


-23.404 


-23.408 




0.60 


1.4 


xl0~ 9 


5.4 


xlO- 12 


0.994 


-23.509 


-23.504 


-23.506 


-23.507 


0.65 


4.2 


xlO" 10 


2.7 


xlO- 10 


0.995 


-23.629 


-23.622 


-23.626 


-23.626 


0.70 


5.6 


xlO- 10 


4.9 


xlO" 11 


0.991 


-23.780 


-23.762 


-23.776 




0.75 


8.2 


xlO" 10 


2.2 


xlO- 10 


0.978 


-23.982 


-23.937 


-23.967 


-23.969 


0.80 


8.1 


xlO" 10 


3.8 


xlO- 10 


0.971 


-24.223 


-24.150 


-24.205 


-24.209 



Table VI. Results from the fit and from the Lanczos-like evaluation for S=l data in Table IV. 
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Fig. 1. A schematic view of the Shastry-Sutherland (SS) model on a 8 x 8 lattice. Solid lines represent the inter-dimer 
coupling (J') and other lines do the intra-dimer coupling (J). 




Fig. 2. Distributions of coefficients for the 4 x 4 SS model in the normal representation (the dashed line) and in the 
restructured representation (the solid line). The ordinate is the absolute values of each coefficient <v while the 
abscissa denotes the basis number i' reordered so that |<v | > \cj> \ holds for i' < f . 
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Fig. 3. Results for 5 = energy difference for the 8 x 8 SS model obtained by the Stochastic State Selection (SSS) 
method with the fit and the Operator Variational (OV) method. 17 ^ Here Eo = — 24. 0J is the ground state energy 
in the singlet dimer phase. 
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Fig. 4. Results for S = 1 energy difference for the 8 x 8 SS model obtained by the SSS method with the fit and the 
OV method. 17 - ) The dotted line is the spin gap estimated by the perturbation theory up to the fifth order. 10 - 1 



